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1 Introduction 



Recent numerical and analytical work on weak first-order phase transitions has shown 
that there may be significant changes to the standard theory of phase transitions as a 
result of thermal fiuctuations Instead of a smooth homogeneous background 

there may be a significant amount of phase mixing due to the existence of subcritical 
bubbles. Though not without controversy these findings suggest we may need 

to reevaluate the standard theory of nucleation of critical bubbles because in very 
weak first-order phase transitions the standard assumption of only small amplitude 
fluctuations breaks down. As the phase transition increases in strength we expect 
the role of fluctuations to diminish until the approximations made for homogeneous 
nucleation become applicable. The effects of thermal fluctuations may extend beyond 
the regime of nucleation however, and alter the dynamics of bubbles as they expand. 
Due to the complex nature of the system, analytic investigations of this question would 
be difficult. On the other hand, numerical simulations in 3-1-1 dimensions would be 
very computationally extensive. We therefore address the problem of dynamics by 
undertaking numerical simulations in 2-1-1 dimensions. Our primary motivation is to 
gain at least a qualitative understanding of how thermal fiuctuations may affect the 
electroweak phase transition in the early universe and the consequent ramifications 
on baryogenesis at this epoch. 

Work on the topic of electroweak baryogenesis has, in general, been concentrated 
in three areas: the form of the effective potential [^j, the dynamics of the transition 
P],[|10|-[ll4l, and how to calculate the baryon asymmetry |[T^-[^. These three areas 
effectively form a hierarchical structure where the means by which one computes the 
baryon asymmetry depends on the dynamics of the transition, which in turn depends 
on the form of the effective potential. Though we do not yet know what the true 
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effective potential looks like, it nevertheless makes sense to investigate the other two 
areas by assuming certain generalities. That is, we may choose a theory with a scalar 
order parameter responsible for spontaneous symmetry breaking. If we assume a 
potential which may be approximated by 



where the coefficients are temperature dependent, we can compute quantities such 
as the temperature needed to nucleate bubbles of the new phase, the thickness of 
the bubble walls, and the rate at which the old phase is converted to new phase 
as a function of the coefficients.^ We require a first-order transition because of the 
third of Sakharov's conditions for baryogenesis, the lack of thermal equilibrium. Be- 
cause the cooling rate of the universe is extremely slow at the electroweak transition, 
depending only on the rate of expansion, the cosmological fluid maintains thermal 
equilibrium. However, by allowing a bubble of stable vacuum to appear within the 
metastable vacuum, a first-order transition gives rise to out of equilibrium processes 
in the neighborhood of the bubble wall. It is only within this confined region where 
baryogenesis can take place. 

The order and strength of the electroweak phase transition is the subject of ongo- 
ing research which will only be solved by a calculation of the effective potential which 
is valid to all orders of perturbation theory. In any case we will assume a potential of 
the general form of eq. (|I]). The strength of the transition is determined by a ratio of 
the coefficients in eq. (|l|), the value of which is temperature and model dependent. In 
the minimal standard model for the electroweak theory the quark, gauge boson and 

Higgs masses determine the height of the energy barrier separating the two minima of 

-'^In fact, the rate of phase conversion depends on the velocity of expansion of the bubble walls, 
which depends not only on the strength of the transition but also on the interaction of the wall with 
the cosmological plasma. 
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the effective potential. If the Higgs mass is small the transition is strongly first-order, 
whereas a large Higgs mass results in a transition which is at most weakly first-order 
and possibly second-order. While most work on phase transition dynamics has been 
done for a strong first-order transition, the experimental lower bound on the Higgs 
mass of 65 GeV PT| has all but ruled out this regime in the context of the minimal 



standard model. More recently, there has been interest in the dynamics of very weak 
first-order transitions [0-0. In this case the phase transition may be completed 
through the process of phase mixing, whereby subcritical bubbles effectively restore 
symmetry, rather than the conventional nucleation of critical bubbles. However, it 
is also possible that the dynamics lie in the intermediate region between a strong 
transition — nucleation and expansion of bubbles in a homogeneous background — 
and a weak transition — phase mixing and domain coarsening. 

This work focuses on the dynamics of a transition which lie in this intermediate 
region. In this case fluctuations will not be strong enough to completely dominate 
the transition and restore thermal equilibrium, yet may have an effect on processes 
such as the nucleation and expansion of bubbles. Though we will not address the 
issue of nucleation (see refs. 0-0] for discussions on this topic), we will investigate 
the ramifications of fluctuations on the speed and structure of the bubble wall and 
what the implications are on the standard mechanisms for generating excess baryons 
at the electroweak phase transition. These mechanisms depend strongly on particle 
interactions with the bubble wall, and calculating the resulting baryon abundance 
requires knowledge of the details of the transition, including the wall's speed and 
structure. Some authors have argued that subcritical bubbles are not relevant to 
baryogenesis since fluctuations only have an effect when the phase transition is only 



very weakly first order, too weak to be of interest to baryogenesis [0|,p|,[15|. These 



arguments, however, deal primarily with the question of whether the phase transition 



proceeds through the nucleation of critical bubbles. The issue addressed in here is 
whether the fluctuations can have a significant effect on critical bubbles after they 
nucleate. 

This issue is related to one which has recently become of interest in condensed 
matter physics. Specifically, there has been an effort to understand how noise affects 
the properties of front propagation in a phase transition [^]. The results of these 
studies has varied depending on the specifics of the model in question. Here, though 
fluctuations play the role of the noise in the phase transition, the system we consider 
and the model we use to analyze it is significantly different from those evaluated 
within the context of condensed matter. 

Once one understands the dynamics of the electroweak phase transition it becomes 
possible to calculate the baryon asymmetry. Traditionally, baryogenesis has been 
investigated in two different limits, thin walls or thick walls. In the thin wall regime 
the wall is thin compared to the mean free paths of the particles, which behave as if 
they were scattering off a potential barrier with CP violating reflection coefficients. 
The reflected charge results in a baryon asymmetry in the region preceding the phase 
boundary. In the thick wall case the plasma is treated as though it were in quasistatic 
thermal equilibrium. Chemical potentials are introduced for quantities which vary 
slower than the time it takes for the bubble wall to pass, and baryogenesis is the 
result of a change in the CP-violating phase, which acts like a chemical potential 



for baryon number |T^. In both cases the transition from one vacuum state to the 
other is smooth, with the order parameter varying monotonically. If fluctuations play 
a significant role it is likely that there would be deviations away from the smooth 
transition model. It is unclear, though, for what transition strength these deviations 
can be ignored. Furthermore, our assumptions about how thick the wall is for a given 
transition strength may have to be modified. 



The plan of this article is as follows. In Section 2 we discuss the potential we use in 
our simulations and its temperature dependence. Section 3 contains the evaluation of 
the nucleation temperature and expansion velocity of critical bubbles in the standard 
homogeneous background model in 2+1 dimensions. We elaborate on our model for 
the phase transition with thermal fluctuations in Section 4, covering such issues as 
the equation of motion, lattice considerations, and specific results. Finally, in the 
conclusion, we discuss ramifications to electroweak baryogenesis as well as possible 
improvements for future work. 

2 The Standard Scenario 
2.1 The Potential 

The potential we select deliberately resembles the temperature dependent electroweak 
effective potential for the minimal standard model. 



The parameters a and A determine the strength of the phase transition and in the 
minimal standard model are related to the gauge boson masses and the Higgs mass 
respectively. The application of this potential, however, may be more general and 
useful in studies of first-order transitions. Due to the 2+1 dimensional nature of the 
simulations,!^ has a mass dimension of [M^] with being [M^/^]. Furthermore, we 
do not associate particular values of a and A with particle masses, but rather con- 
centrate on which regions of parameter space constitute weak and strong transitions 
as determined in Ref. 0. 

We will find it useful to change from dimensional to dimensionless variables as 
follows. 




(2) 



X ^ x/ \faT2 



(3) 
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t t/v^Ts (4) 

- XTI" (5) 

T ^ dT2. (6) 
The potential is now, 

V{x, 0) = \(0'- - \^0X' + ^Ax^ (7) 

where a = aa-^T^^''^ and A = Aa^^Ta^^ For ^ < 0^ = (1 - aV4A)-i/2 ^j^g potential 

has one maximum X- and two minima Xo? X+ located at Xo = and x± = a^/2A(l ± 



1 - 4A(1 - 6'-2)/a2). At the critical temperature 6*^ = (1 - 2aV9A)"^/^ the two 
minima are degenerate with the minimum at x+ being the global minimum for 6 < 6c. 
At 9 = 1 the minimum at Xo disappears. Figure 1 shows the potential for A = 0.1, 
a = 0.4, 9 = 9^- 0.0025 = 1.2432. 

2.2 Dynamics 

The dynamics of the phase transition depends to a large degree on the amount of 
supercooling undergone before the nucleation of bubbles. If the supercooling is large, 
as in the case of a strong first-order transition, the bubbles expand very rapidly and 
have a relatively thick wall. In a weak transition the free energy of a bubble is 
minimized for a thin wall and the small cooling results in a slowly moving wall. The 
first step, then, is to determine the temperature at which nucleation takes place. 

As the universe cools below the critical temperature the symmetric vacuum be- 
comes metastable with a finite probability of decaying into the asymmetric stable 
vacuum. The theory of bubble nucleation from a metastable to a stable state was 
developed by Langer and later applied to cosmological phase transitions ||2^,[l5 . 
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The equation describing the rate of nucleation per volume in 3+1 dimensions, how- 
ever, must be changed to account for the different scahng in a 2+1-dimensional modeL 
Recall that in 2-dimensions volume scales as and not T^^. The nucleation rate 
per volume is 

r/V = AT^e-^^/^, (8) 

where Fc is the free energy of a critical bubble. The rate is dominated by the expo- 
nential; hence, the exact value of A is not very important and we set it equal to one. 
The nucleation temperature is given by the temperature at which the probability of 
nucleating a critical bubble inside a horizon approaches one. We must, therefore, 
determine the size of the horizon as a function of temperature. In a 2+1 dimensional 
universe during the radiation dominated era, the energy density of the universe p 
varies with the scale factor a according to p oc a~^. The time-temperature relation is 
then 

where mp\ = 1.22 x lO^^GeV is the Planck mass and ^ ~ 1/30 ^8[. The volume 
inside the horizon at a temperature T is 

Vh « 4^^^. (10) 

From this we can write down the probability of nucleating a bubble inside a causal 
volume at a temperature T, 

dP ^ VVndt ^ ee (^)'^'e-^=/^f^. (11) 

We define the nucleation temperature as the temperature for which the total proba- 
bility of having nucleated a bubble approaches one, 

l-4e^(fi)'''e-^^/^", (12) 



where T„ is the nucleation temperature. Estimating T„ is made much simpler by 
approximating T„ in the prefactor of the right hand side of Eq. (0) as the critical 
temperature. Because of the exponential, the final answer is not very sensitive to the 
value of the critical temperature chosen. A value of = 100 GeV yields 

53 ^ FjTn. (13) 

In order to calculate the free energy of a vacuum bubble, we choose the energy of 
the meta-stable vacuum as zero, V^((0) = 0) = 0. The excess free energy of a bubble 
is 

The first term in the integral represents the surface energy of the bubble, while the 
second term is the volume energy, coming from the difference in free energy inside 
and outside the bubble. The free energy, in general, must be found numerically. The 
difficulty in calculating F arises because one needs to know the value of the field at 
all points in space, which in general may not take a simple functional form. Under 
certain limits, however, approximations prove fairly accurate. In the case of extremely 
small supercooling the wall approaches the well known kink solution, the "thin-wall" 
approximation. As the temperature drops further below the critical temperature 
the wall shape is reasonably approximated by a Gaussian ||10|. At temperatures 
appropriate to nucleation, neither approximation turns out to be particularly good. 
In spite of the fact that the thin-wall approximation fails to accurately predict the 
free energy associated with the nucleation of critical bubbles, we use it here as a way 
of estimating a value for T„. One should bear in mind, however, that the free energy 
of the true critical bubble solution is larger than the free energy in the thin-wall case 
at a given temperature. As a result the nucleation temperature in the thin-wall case 
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is higher, i.e., less supercoohng, than in the true case. We, however, are not most 
interested in what the actual nucleation temperature is, but rather in the dynamics of 
bubble expansion. Within this context the thin-wall approximation gives a reasonable 
estimate. Also, recent studies have shown that the actual nucleation temperature may 
be significantly higher than what one obtains by the standard method [^], |0. The 
basic idea is that in a weak first-order transition large amplitude fiuctuations cause 
the energy density of the metastable vacuum to shift away from V{0). Instead one 
must include a non-perturbative correction to the free energy of a critical bubble 
which has the effect of raising the nucleation temperature. 

To calculate the free energy in the thin- wall limit, when supercooling is small, 
we use a perturbative expansion in A = 1 — 6/6c, the amount of cooling below the 
critical temperature. To first order the potential is 



A bubble of true vacuum which is just large enough to grow satisfies the static solution 
to the equations of motion. 



dr"^ r dr dx 

In the thin-wall limit the spatial derivatives are small except in the bubble wall. 
Furthermore, when A is small the radius of a critical bubble becomes very large and 
the first order derivative term in Eq. (|16D becomes negligible. This gives 




(15) 



X Idx dV 




dr 



(17) 



The free energy of a bubble of radius R is 
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where x+ represents the value of the field of the true vacuum and is given by 

As a function of the radius R and the amount of supercooling A the free energy is 

F(fl,A)«— y2A^^j - — A. (20) 

By taking the derivative with respect to R we get the radius of a critical bubble 
i?c = a;6'c/(18v^2AA). Since A is small, the nucleation temperature is approximately 
the critical temperature, and after equating with Eq.(|l3|) we get for the critical free 
energy to temperature ratio 



Qn 1458 A^A 



53 (21) 



or, 



A ^ (4.07 X 10-^)^. (22) 

Once a critical bubble is nucleated it begins to expand because the free energy lost 
due to the interior being in the lower energy vacuum offsets the gain in surface energy. 
The wall quickly accelerates, and in the case of a vacuum transition, approaches the 
speed of light within a few wall widths pS. In the more realistic case where the 



wall is expanding through a plasma, the wall experiences an opposing force due to 
interactions with particles and reaches a terminal velocity. Calculating the terminal 
velocity proves to be a difficult problem because it depends on detailed interactions of 
the Higgs field with the plasma. Furthermore, the size of the damping also depends on 
the thickness of the wall relative to the mean free paths of the particles in the plasma 
T0| . Simply, this results because a thin wall causes particles to make the transition 



into the true vacuum state quicker than the time it takes for them to equilibrate. A 
thicker wall allows the particles to maintain quasistatic thermal equilibrium, because 

10 



the thermalization rate is faster than the rate at which the Higgs field changes, but not 
chemical equilibrium, because some particle interaction rates occur slowly resulting 
in departures from equilibrium populations. The result is a different damping force 
depending on the regime. In any case, we assume here that the damping force can 
be modeled by a velocity dependent term in the equation of motion of the Higgs 
field, where the magnitude of the proportionality constant determines the terminal 
velocity, 

(fx , dx „2 dV 

^ + ^^-^^ = -^- ^^'^ 
In the frame moving with the wall we can change coordinates to r = 7(0; — vt), 
simplifying to the case of a very large bubble so that we may treat the problem as 
effectively one dimensional. The equation of motion then takes the form 

d\ , dx dV 

The boundary conditions on x state that the derivative must vanish far from the wall. 
We can then integrate to obtain 

"'^'^ IZ (^) = " ^^^^ 

= -Vix+) (26) 

Within the validity of the thin-wall approximation, we may perform the integral 
analytically to obtain 

Iviv^xl = Vix+)- (27) 

Several authors |§| , , [ll]] , [Q have calculated the velocity of the bubble wall in both 
the thick and thin case. The results for a thin wall are f ~ 0.1 and v ~ 0.2-0.6 for a 
thick wall. Given a particular set of parameters for the potential, then, we can find 
an appropriate value for rj. Here we say that rj is generally (9(0.1). 
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Up to this point our treatment of the bubble has been from a purely field theoretic 
point of view. We describe the dynamics of the Higgs field as a scalar field in a tem- 
perature dependent potential. Particles in the plasma scatter off the field effectively 
creating a damping force. A more macroscopic view of the phase transition describes 
the plasma through hydrodynamics and the bubble wall as an effective combustion 
front. In this case there are several effects which arise which affect the dynamics 
of the bubble expansion which are not otherwise evident in the field theoretic point 
of view. The velocity of the wall, for instance, depends not only on the damping 
coefficient rj, but also on the ability of the fluid to conduct heat away from the wall 
and the resulting small temperature deviations created through the release of latent 



heat [|TT],|jTj]. Furthermore, fluid dynamical instabilities can arise both in the plasma 
and the bubble wall. Perturbations in the wall may grow exponentially depending on 
the size of the perturbation and the strength of the front [^-[Q. The net effect is 
that bubbles might not grow in a spherically symmetric way. Nevertheless, including 
these effects is beyond the scope of this work and we consider here only the simplified 
scenario of a scalar field in a potential well with damping. 

3 Phase Transitions with Fluctuations 

Thermal fluctuations play an essential role in the phase transition, so it is instructive 
to discuss briefly how they enter into the physical picture. Without fluctuations the 
field remains trapped in the metastable vacuum until the potential barrier vanishes 
at which time the field rolls down the potential into the true vacuum. In a thermal 
environment, though, it is only the expectation value of the field which is in the 
metastable state. There is a finite probability that the field can take on a value beyond 
the potential barrier and thus make the transition into the true vacuum. Hence, the 
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nucleation and expansion of bubbles. This method, however, traces the evolution of 
the field using the equations of motion at zero temperature often with a damping term 
to simulate the dissipation due to the thermal environment, with finite temperature 
effects limited to corrections to the effective potential. This approach is justifiable 
if one wishes to consider only the behavior of the expectation value of the field, 
where one has implicitly assumed that the variance of the field value is very small. 
This assumption should hold if the minima in the effective potential are sufficiently 
deep that thermal fluctuations away from the minima are strongly suppressed. Even 
in this case, however, the dynamics of the field may not be adequately described 
by the equations of motion. For instance, if the field is initially located at a local 
minimum of the potential, which then becomes a local maximum in such a way 
that the slope remains zero, the field remains in this unstable extremum indefinitely. 
Realistically, what happens is that fluctuations dislodge the field from the extremum, 
which subsequently evolves according to the equations of motion. 

The question of how one is to model the thermal fiuctuations in field theory 
nevertheless remains. In classical statistical mechanics one treats the problem of 
thermal fiuctuations by identifying a system and a heat bath and choosing a model 
for the couphng between the two and solving for the behavior of the system. Finding a 
solution is greatly simplified when the reservoir obeys particular criteria. Specifically, 
if the reservoir has infinite specific heat and a relaxation time much shorter than 
that of the system, it remains in thermal equilibrium even though it interacts with 
a system which may not be. When these conditions are satisfied one may ignore the 
dynamics of the reservoir in favor of the dynamics of the system. This paradigm 
provides a method for calculating quantities such as equilibration time scales of a 
system and transition probabilities when the system is in equilibrium. In the case of 
Brownian motion, for example, where the system is a small macroscopic particle and 
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the bath is the fluid in which it rests, the Langevin equation is a natural outcome of 
the fluctuating thermal force of the bath on the system. 

In the case of field theory, however, the dichotomy of the system and bath is 
less clear. If there is only one self-interacting field one can decompose the dynamics 
into short and long-wavelength modes which operate on different time scales. The 
short wavelength modes respond much more quickly and can serve as the thermal 
bath while we define the system as those modes whose wavelengths are larger than 
some critical value; non-linear interactions couple the system to the bath. Another 
approach is to have a second field which acts as the bath to the first field. In either 
case the form of the coupling determines the nature of the fluctuations which the 
system experiences. Below, we model the fluctuations as white noise (uncorrelated) 
by adding a stochastic term to the equation of motion. In general, the noise which 
the system experiences may be significantly more complicated, as is the case in some 
simple cases which have been studied ||3^,[^. We justify our choice by noting that 
it is not yet clear how these findings might change given the relatively complicated 
environment of the early universe. Furthermore, in at least on study |Q the authors 
found that in the high temperature limit the noise becomes white. 

The Langevin equation is a popular way to model phenomenologically a system 
with thermal fluctuations, though, as mentioned above, the actual dynamics of the 
noise may be different from the simple model of additive white noise. Nevertheless, 
the crucial question is deciding how one is going to model the essential physics of 
the system. Certainly, if we knew what the true equations of motion for the Higgs 
and particle fields were in the early universe, we could use those equations. This, 
however, is not the case. At the time of the electroweak phase transition, aside from 
the Higgs, there were many other particles, and writing the complete dynamics is, if 
not impossible, extremely difficult. Given this situation, using the Langevin equation 
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with additive white noise is a reasonable place to start in the investigation of how 
fluctuations may affect the dynamics of the transition, and we use it here. Also, 
numerical studies of nucleation of critical bubbles in 1+1 and 2+1 dimensions carried 
out using the Langevin equation to model thermal fluctuations have demonstrated 
good agreement with classical nucleation theory [^],|^ . 



3.1 The Model 

As described above the coupling of the field with the thermal bath is modeled by 
a Langevin equation. The equation of motion, then, in terms of the dimensionless 
variables is 

where we have defined C, as the dimensionless stochastic noise. The fluctuation- 
dissipation theorem relates the noise to the dissipation coefficient rj by 

(e(x, t)e(x', t')) = 2i^e5{t - t')5{^ - x'). (29) 

In discrete form, ^ is 

^(Xj,t„) = S,i,n (30) 

Qr.,n (31) 



where Qi^n is a unit variance Gaussian random number at each point on the lattice 
and 5t and 5x are the lattice spacing in the time and spatial directions, respectively. 
We integrate the equation of motion forward in time using a second-order leap-frog 
method. 

We carry out the numerical experiment by inputting initial conditions and allowing 
the simulation to run. Since we are interested in the dynamics of the bubble wall 
and not questions of nucleation, we use as initial conditions a wall which stretches 
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across the width of the lattice located at the midway point along its length. Half 
of the initial lattice volume, then, is in the asymmetric phase while half is in the 
symmetric phase. The shape of the wall is chosen so that it conforms to the kink 
profile appropriate for a thin wall. This profile is in fact not the one which minimizes 
the free energy, but it is sufficiently close so that the time it takes for the wall to relax 
to its correct form is much smaller than the run time. By making the wall stretch 
across the width of the lattice we effectively model a large bubble and thus ignore 
the initial expansion stage following nucleation. The wall, however, is not given any 
initial speed, but must accelerate to its terminal velocity. The time for this to happen 
is also very small. 

The initial configuration of the wall is prepared without noise added in. We 
include the fiuctuations in the simulation only through the dynamics of the equation 
of motion. Though such a configuration is certainly unphysical given the assumption 
that thermal fiuctuations exist, it has no effect on the outcome of the numerical 
experiment. In each of the two phases there exists a thermal equilibrium distribution 
of the other phase due to fluctuations Given the strengths of the transitions 

we consider, the time it takes to reach this distribution is small. Only during the 
initial stages of the simulation is the bubble wall likely to be affected by particular 
characteristics of the initial conditions. 

A flrst-order phase transition may be classifled into weak and strong transitions 
depending on the height of the barrier separating the minima. In a weak transition 
considerable phase mixing may exist with the transitions proceeding through domain 
coarsening, while in a very strong transition we expect the theory of homogeneous 
nucleation and bubble expansion to be correct. We wish to explore here the intermedi- 
ate region and thus must choose appropriate values for the parameters which describe 
the potential. Based on the flndings in Ref. we select values for the parameters 
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a and A which explore transitions in the strong regime, where a strong transition is 
one where at 6c the symmetric phase comprises at least about 60% of the total area. 
We, therefore, set A = 0.1 and allow a to take on the values 0.4 and 0.5, well into the 
strong regime. 

The potential in Eq. has a direct connection with the electroweak phase tran- 
sition. In the temperature one-loop effective potential a is related to the gauge-boson 
masses while A is determined by the Higgs mass. The values we choose here for these 
parameters, however, should not be construed as exploration of the parameter space 
of masses of these particles. The naive interpretation that the chosen values of a and 
A correspond to masses in the minimal standard model is false in this case because 
the simulations describe dynamics in a 2+1 dimensional world. What constitutes a 
weak and strong transition in 2+1 dimensions is different from the 3+1 dimensional 
case. We choose, therefore, not to make any claims about particle masses, but rather 
focus qualitatively on the effects of fluctuations on the phase transition. 

3.2 The Lattice 

When taking a numerical approach to this problem one must make a choice regard- 
ing the coarse-graining scale as realized through the lattice spacing. Modes with 
wavelengths shorter than the lattice spacing couple to the larger wavelength modes 
through the noise term in the equation of motion. The results from placing the field 
theory on a lattice, then, only apply to the long wavelength modes. When probing 
physics at shorter wavelengths, one must be careful in taking the continuum limit. To 
do so one must include renormalization counterterms in order to obtain the proper 
continuum theory. These issues are discussed in more detail in Ref. |3^. Following 



the prescription in Ref. we set the lattice spacing 6x = 1, which is approximately 
equal to the mean-field correlation length given by = V"{xo, 6c)- 
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In addition to choosing a grid size on the lattice, one must choose the size of a 
time step, the size of the lattice, and the boundary conditions. Ideally a large time 
step would allow one to integrate the equations of motion for a longer period of time; 
however, stability considerations limit the size of 6t. To find an appropriately sized 
step we compared simulation results as a function of different values of St, choosing a 
value in the regime where the results become independent of the magnitude of step. 
In all simulations wc used a value of 6t = 0.2 

As with the choice for St, selecting a lattice size is another exercise in compromises. 
The physics one is attempting to simulate takes place in an effectively infinite volume, 
but one is limited to not just finite lattices but also ones which are fairly small because 
of the need for reasonable integration times. There are dangers, however, in having 
too small a lattice. Strictly speaking, in the context of phase transitions symmetry 
breaking only occurs in the infinite volume limit. Within any finite volume there exists 
the possibility that a fiuctuation will restore the symmetry regardless of the dynamics 
internal to the volume. Since it is precisely these dynamics with which we concern 
ourselves, it is paramount to choose a volume large enough that the probability of 
such an occurrence becomes negligible. Essentially, what happens is that fiuctuations 
in the broken phase may drive the system back to the symmetric phase, even though 
it is energetically unfavorable. If the total volume is large, such fluctuations result 
in only a small region of the total volume having its symmetry restored. The chance 
that this could happen with a large volume is exponentially suppressed because of 
the large amount of energy necessary. We can estimate how large a volume we need 
by calculating the rate for a large fluctuation, 

r - e-^/^ (32) 

The free energy F is given by the change in effective potential energy between the 
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two minima multiplied by the amount of volume in the asymmetric phase. Ignoring 
the gradient portion of the free energy, which only increases F and makes a symmetry 
restoring fluctuation less probable, we have 

J " (^^) 

where V is the volume in the symmetric phase. The probability of having a large 
fluctuation given the size of our lattice and run time is approximately 10~^°. 

Since we are interested in the dynamics of the bubble wall as it expands, we 
also want to make sure that it is unlikely that a critical bubble will nucleate in the 
asymmetric phase within the lattice during a run time. Because the volume under 
consideration is much less that a horizon volume and the time is much less than a 
Hubble time, the probability of this happening is similarly suppressed. The lattice 
size we use in the simulations is — 100, Ly — 50, where is the length in the x 
direction and Ly is the length in the y direction. 

Another consideration lies in the fact that we are modeling an unbounded system 
by a flnite volume with a boundary. Though, as discussed above, we do not expect 
finite size effects to be important, there is a distinct surface in the simulations which is 
unphysical. We circumvent this issue by using periodic boundary conditions in the y 
direction and 'open' boundary conditions in the x direction. The periodic conditions 
allow us to model an essentially planar wall, appropriate for a large bubble. The 
danger, however, with periodic conditions is that long range correlations may be 
induced if the simulation is run for longer than Ly/2, the light crossing time. It turns 
out that this is not a concern here because the presence of dissipation and noise have 
the effect of damping out and swamping any long range 'communication' which might 
otherwise exist. The open boundary conditions consist of assuming that for points 
immediately outside the lattice the field takes on a value equal to the field on the 
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boundary. Though not reahstic, these boundary conditions provide an approximation 
to the unbounded system. Any spurious effects caused by these conditions do not 
extend into the lattice because of dissipation and noise. 

4 Results 

The focus of the numerical experiment is to understand the effect of the fluctuations on 
the rate at which the old phase is converted into the new phase and how that transition 
is made. We expect that for a relatively stronger transition the fluctuations will play 
a fairly minor role, while for a relatively weak transition there may be a marked 
difference from the homogeneous background case. We investigate the rate of phase 
conversion for different transition strengths by allowing a to vary while holding all 
other parameters constant and comparing the results to the case where there are no 
fluctuations. The values of a we choose (0.4 and 0.5) place the transition in the strong 
regime as defined in Ref. 0. In that study the author found that the change from 
a weak to a strong transition is itself a second-order phase transition with a critical 
value of = 0.36. We can characterize the strength of the transition by /+, the 
fraction of volume which fluctuates from the symmetric phase beyond the maximum 
of the potential barrier. At the critical value this fraction is 42%, while for a = 0.4, 
/+ ~ 6% and for a = 0.5, /+ ~ 0.1%. Furthermore, whereas in the homogeneous 
background case the field in a region of space smoothly makes the transition from the 
symmetric to asymmetric phase, we do not expect this to happen when the amplitude 
of the noise becomes large. 

In order to calculate the rate of phase transformation for a particular transition 
we first introduce a definition which will allow us to determine which phase the field 
is in. We label the field at a particular lattice point as being in either the 0-phase if 
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X < X_ (i.e. to the left of the maximum) or the + phase ii x > X- (i-G- to the right 
of the maximum). This allows us to determine what fraction of the total volume is 
in each of the two phases. Because we place the bubble wall through the center of 
the lattice approximately half of the initial volume is in each phase with deviations 
away from half due to fluctuations. Lattice points move on average from the 0-phase 
to the + phase as the bubble expands. In Figures 2-3 we plot the position of the wall 
as a function of time for a phase transition with and without fluctuations for a = 0.4 
and a — 0.5 with rj — 0.2. The wall position for the case with fluctuations as shown 
by the short-dashed line is an ensemble average of 200 separate trials while the long- 
dashed lines show the one standard deviation width of the distribution. The dotted 
line shows one realization. The solid line is the position for a phase transition without 
fluctuations, where the step-like nature is due to the discreteness of the lattice. For 
the case with fluctuations the position does not mean that the midpoint of the wall is 
located at that particular x-position all the way across the lattice rather it represents 
an average position at that point in time. In fact there are regions both in front 
of and behind the wall which belong to the opposite phase resulting in a somewhat 
amorphous boundary. Figs. 4a-d show contour plots of the fleld on the lattice for 
a — 0.4, 0.5 with rj — 0.1, 0.5 with contours at x_/2 (dotted), X- (solid), 3x-/2 
(dashed). We see signiflcant distortions in the wall away from planar in the weaker 
transitions whereas the stronger transitions approach the plane solution, though other 
structure is still evident. The rate at which phase is converted from the to the -|- 
state is also clearly elevated relative to the transformation without fluctuations. In 
Table 1 we show the rate of phase conversion, essentially the velocity of the wall, for 
the different transition parameters. The right-hand columns for a given value of a and 
rj show the average wall velocity for the fluctuating background case along with the 
one sigma width; the left-hand columns show the wall velocity for the homogeneous 
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background case. In the limit of a very strong transition we can see that the rate of 
phase conversion approaches the homogeneous theory solution. 

Also particularly relevant to electroweak baryogenesis is the actual structure of 
the wall. Whether and how many baryons are created depends on how particles, 
which are interacting with the changing Higgs field, make the transition from the 
to the + state. Though this depends on quantities such as the velocity of the particle 
relative to the wall, which we don't calculate here, we can still get an idea of what 
a particle "sees" as it moves from one state to the other by taking cross sections 
of the wall. In Figs. 5-6 we plot the values of the field x as a function of x at a 
particular time for a = 0.4 and a = 0.5 with 1] = 0.1 for both the homogeneous 
and fluctuating background theories. While the wall thickness in the homogeneous 
theory can be clearly defined, it is less apparent in the fluctuating case. We see that 
averaged over 200 simulations the wall is thicker relative to the case where there are 
no fluctuations, though the increase is only moderate for the transition strengths used 
here and approaches the homogeneous background solution for a stronger transition. 

Figures 7-8 follow the behavior of a lattice point as a function of time. The 
boundaries of when the transition from the to + state begin and end is less clear in 
the fluctuating case as compared to the homogeneous case. The results from averaging 
over many simulations indicates that in general a lattice point makes the transition 
more slowly than if there were no fluctuations. This results because in the vicinity 
of the bubble wall a lattice point may undergo several "transitions" before finally 
reaching equilibrium in the + phase. 
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5 Conclusions 



We see that quantities such as the rate of phase conversion and the structure of the 
bubble wall can be different for phase transitions with fluctuations as compared to 
the standard homogeneous background model. The speed of the bubble wall may 
be increased by a factor of about two for a "mildly" strong transition and probably 
even more for weaker transitions. Physically, this is a reasonable expectation; the 
wall effectively "swallows" the field fluctuations, moving forward as it does so. In 
a transition with more large amplitude fluctuations this swallowing effect is more 
pronounced. Though the rate at which phase is converted from the to the + phase 
is increased globally, locally the transition is not necessarily well defined and more 
gradual. Furthermore, we see that this effect is noticeable even though only a small 
percentage (about 6%) of the symmetric phase may fluctuate beyond the barrier at 
any given time. As the amplitude of the fluctuations decrease the bubble behaves as 
predicted by the homogeneous background theory. That there are large amplitude 
fluctuations which can roughly "mask" the transition is not surprising; in fact, this is 
entirely determined by the stochastic noise term in the equation of motion. What is 
interesting is that the noise should have such a significant effect on a transition which 
one would consider "strong" by the definition given in Ref. 

It is not necessarily evident that random fluctuations should result in an increase 
in the velocity of the bubble wall. Though fluctuations from the to + phase might 
speed up the wall, one might also expect that fluctuations from the + to phase 
would have the opposite effect, resulting in no net change. It is the asymmetry in the 
potential, however, which prevents this from happening. Because fluctuations are a 
result of a random impulse on the field at a particular point, a fluctuation which drives 
the field from the to + phase is more likely to approach or exceed the maximum of 
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the potential than those that go from the + to phase. The dominant effect on the 
dynamics, then, is to enhance to to + transition. 

The imphcations for baryogenesis stem from the fact that models rely on particular 
assumptions of the bubble wall structure, e.g., a smoothly varying monotonic order 



parameter. The standard picture |T3] states that, depending on the thickness of the 
wall, one selects a mechanism with which to calculate the baryon number generated. 
In order for this to work, however, it is necessary that the homogeneous background 
theory of bubble nucleation and expansion be valid. Previous studies have called into 
question this paradigm by pointing out that subcritical bubbles may affect the initial 
conditions of the transition as well as bubble nucleation in a weak first-order 

phase transition. Here we investigate transitions which are stronger than those in 
0, but still within the regime where nucleation may be affected by the presence of 
subcritical bubbles B. Although the electroweak transition is most likely weak in 



the minimal standard model, authors |0, P, ||I5[ have argued that baryogenesis is 
likely only for a sufficiently strong transition where the effects of phase mixing due 
to thermal fluctuations would be negligible. What we have shown here is that the 
realm of the phase transition for which fluctuations may play a significant role is 
larger than what was expected. Their effect is not limited to possible alterations in 
the picture of homogeneous nucleation or degree of phase mixing, but also includes 
the dynamics subsequent to nucleation. For instance, not only are the wall thickness 
and velocity larger, but the path from the symmetric phase to the asymmetric phase 
is hardly smooth or monotonic. Any model for electroweak baryogenesis will clearly 
have to take the stochastic nature of the dynamics into account when considering the 
interaction of particles with the Higgs field. The results here are not meant to be 
quantitative, but they do demonstrate the importance of fluctuations in the dynamics 
of first-order phase transitions. 
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Finally, we point out some of the limitations of this work. Ultimately, we would 
like to be able to obtain a complete picture of the electroweak phase transition in 
order to determine the viability of generating the baryon asymmetry of the universe 
at this scale. What we have done here is to examine the role thermal fluctuations are 
likely to have on the transition dynamics, and what the implications are on models 
of baryogenesis. However, this study was limited to 2+1 dimensions rather than the 
3+1 dimensions of the real world. Extending this work to higher dimensions would be 
a valuable contribution to our understanding of this problem. Also, we have assumed 
here that one can model the thermal fluctuations through stochastic white noise, as 
one does when studying phenomena such as Brownian motion. Realistically, this 
is only an approximation of more complicated couplings between system and bath. 
Indeed, how one is to divide the physical environment into system and bath and the 
nature of the noise that results is a problem currently under active investigation ||32|| , 
p3| . Lastly, the dynamics of the cosmological fluid have been omitted. Not only 
should the fluctuations have an effect on the fluid, but, as already noted, the fluid 
itself plays a role in the overall dynamics. A more complete treatment would include 
fluid velocity, pressure and temperature as well as parameters for heat capacity and 
conductivity in the simulations. In spite of these caveats, we believe that our flndings 
are indicative of what we may eventually flnd to be the "true" dynamics. 
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Figure Captions 

Figure 1. Potential for A = 0.1, a = 0.4, 9^9^- 0.0025 = 1.2432. 

Figure 2. Average wall position as a function of time for a — 0.4, 77 = 0.2. Shown 
are a phase transition without fluctuations (solid line), with fluctuations (dotted 
line), average over 200 simulations (short-dashed line), and one-sigma width 
(long-dashed line). 

Figure 3. As in Fig. 2 with a — 0.5. 

Figure 4. Contour plot of field x the lattice with contours x_/2 (dotted), X- 
(solid), 3x-/2 (dashed), a) a = 0.4, 7] = 0.1. b) o; = 0.4, 7] = 0.5. c) o; = 0.5, 
77-0.1. d) a = 0.5, 77 = 0.1. 

Figure 5. Wall profile for a = 0.4, rj = 0.1. Shown are simulation resluts for a 
phase transition without fluctuations (solid line) , with fluctuations (dotted line) , 
average over 200 simulations (short-dashed line), and one-sigma width (long- 
dashed line). 

Figure 6. As in Fig. 5 with a — 0.5. 

Figure 7. Lattice point history for a — 0.4, 77 = 0.1. Shown are simulation results 
for a phase transition without fluctuations (solid line) , with fluctuations (dotted 
line), average over 200 simulations (short-dashed line), and one-sigma width 
(long-dashed line). 

Figure 8. As in Fig. 7 with a — 0.5. 
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Table Caption 

Table 1. Wall velocities for varying parameter values where Vh represents the velocity 
in simulations without fluctuations and Vf represents the velocity in simulations with 
fluctuations. 
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7/ = 0.1 


77 = 0.2 


77 = 0.5 
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Vh 


^/ 


Vh 




Vh 




0.4 


0.15 


0.24±0.06 


0.08 


0.16±0.05 


0.035 


0.077±0.036 


0.5 


0.4 


0.40±0.06 


0.21 


0.25±0.02 


0.09 


0.12±0.01 
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